The Phase Diagram of Two Color QCD 
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Abstract. I present recent results from lattice simulations of SU(2) gauge theory with Nf = 2 
Wilson quark flavors, at non-zero quark chemical potential fi. The thermodynamic equation of 
state is discussed along with the nature of the high density matter which forms. It is conjectured 
that deconfmement may mean different things for bulk and Fermi surface phenomena. 



1. Why Two Colors? 

In this talk I will give an update on our project to simulate SU(2) gauge theory with Nf = 2 
quark flavors on cold lattices with quark chemical potential \x ^ [H02G]. Because of the Sign 
Problem, QC2D offers currently the best prospect of using lattice simulations to study gauge 
theories at non-zero baryon charge density. In particular, we can explore the systematics of 
the lattice approach in this unexplored physical regime. There is good news: unlike the case of 
hot QCD, it is possible to perform a scan in /x at fixed cutoff, and the primary thermodynamic 
observable, the quark number density n q = —df/dfi, as a component of a conserved current 
has no quantum corrections. The bad news is that both UV and IR artifacts are significant and 
complicated; only recently have we begun to make progress in disentangling them [3]. 

A second, more theoretical motivation is that QC2D with /i ^ offers the chance to explore 
color deconfinement in a new physical regime, complementary to the more usually studied 
transition at T > 0. There are interesting differences which will challenge us to refine the 
language used to discuss this essential feature of non-abelian gauge theories. 



2. The Simulation 

We define our lattice theory with the Wilson plaquette action and unimproved Wilson fermions 
with hopping parameter k. Chemical potential is introduced via the standard prescription of 
weighting forward/backward temporal hops by e ±fJ,a . Reality of the fermion measure follows 
from detM = detT2M*T2 where the Pauli matrix T2 acts on color. Our only innovation [1] 
is the introduction of a scalar, isoscalar and gauge-invariant diquark source term jqq = 
jK(ip2~Cj5T2tpi — '0iC'75T2V>2 T )> where subscripts denote flavor. Setting j / mitigates the 
IR fluctuations due to Goldstone modes associated with superfluidity at non-zero quark density. 
Indeed, j is nothing but a Majorana mass for the quarks. Fig. [1] shows how the computational 
effort, measured by both the number of congrad iterations and the HMC timestep dt required 
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Figure 1. Computational effort needed for 
QC 2 D. 
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Figure 2. n q /n^ B and pjp for ja = 0.04 



on 12 d x 24. 



for decent acceptance varies with /i and j. Despite the absence of a Sign Problem, studies of 
dense matter still require an order of magnitude more cpu than that needed for vacuum QCD. 

By now we have accumulated matched ensembles on a coarse 8 3 x 16 lattice with a = 
0.229(3)fm p] and a finer 12 3 x 24 lattice with a = 0.178(6)fm [2], with scale set by assuming 
the two color string tension is (440MeV) 2 . Both lattices have temperature T ~ 50MeV, making 
them "cold" in QCD terms. The quark mass is relatively heavy, corresponding to m n /m p ~ 0.8 
- up to now we have not considered the chiral limit to be so important, because effects due 
tom/0 could naively be expected to lie deep at the bottom of the Fermi Sea. In our most 
recent simulations [3] we have explored the systematics both of varying the diquark source 
ja = 0.04, . . . , 0.02 and varying T, with N T = 24, . . . , 8 corresponding to T = 47, . . . , 141MeV. 



3. Equation of State 

Fig. [2] plots n q /rig B as a function of p,a, calculated on the fine lattice with ja = 0.04 [2j. Here 
n q B denotes the "Stefan-Boltzmann" (SB) result evaluated for free Wilson fermions on the same 
volume; taking the ratio is a first step towards correcting for discretisation and finite volume 
artifacts. These are considerable, as shown by the departures from unity in Fig. [5j For [ia 
greater than unity, the dominant error is due to UV artifacts, but for smaller u the disparity 
between the two curves points to the IR; indeed, the oscillatory behaviour is due to departures 
of the Fermi surface from sphericity as T — >■ due to the discretisation of momentum space [1] . 
Fig. [3] illustrates the difficulty of disentangling these two effects, updating Fig. [2] for a range 
of temperatures with the limit j — > taken. Fig. [J] presents the same data normalised by 
the continuum SB result NfN c (fj, 3 /n 2 + / uT 2 )/3. The apparent peak at \ia « 0.4, previously 
identified as a signal of Bose-Einstein condensation (BEC) of tightly bound diquarks [HE], is 
most likely an IR artifact due to the dip in the dashed red free fermion curve of Fig. [5j On the 
other hand the difference in vertical scales shows that UV artifacts still need to be corrected. 

Another previously neglected factor is the j ' — > limit, which appears to be a much more 
important effect for interacting quarks; setting enhances their tendency to form pairs, in 

effect promoting BEC formation which could deform the results. Taking all factors into account, 
we conclude the approximate plateau in the range 0.4 < /ia < 0.7 is consistent with n q /n^ B ~ 1, 
rather than 1.4 as suggested by Fig. [2] The same gross features are also apparent in a calculation 
of the pressure p = f 1 * n q d\i\ Fig. [6] shows the result of three different approaches to correcting 
for lattice artifacts at T = 47MeV [3]. 

Figs. I3|4l and IB1 support an interpretation in terms of four distinct regimes separated by three 
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Figure 3. n q /n q for j — > using lattice 
free fermions. 



Figure 4. n q /n q B for j — > using continuum 
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Figure 5. 37r 2 n^ s /NfN c /j, 3 vs. //a on two 
different lattice volumes. 
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Figure 6. p/p s on 12 3 x 24 for j -> usme 
three different artifact correction schemes. 



thresholds. The first is the onset threshold, taking place at \x = m 7T /2 ~ 0.32a _1 ~ 360MeV, 
separating the vacuum state from a non-zero quark density n q > 0. In mean field theory this 
is predicted to be a second order phase transition [5], unlike the first order transition expected 
in QCD. The non-monotonic behaviour of n q /n q above onset due to a BEC, predicted in 
mean field-theory, and reported in [HE] is, however, no longer a robust feature of the numerical 
results; it may well be that simulations far closer to the chiral limit are needed to expose this 
behaviour. Next, at hq ~ 0.5a" 1 ~ 530MeV is a crossover to a regime where both n q and p are 
approximately equal to their SB values, suggestive that a degenerate system of quarks with a 

1/3 

well-defined Fermi momentum kp oc n q and Fermi energy Ep ~ fx has formed. For reasons that 
will be elaborated we refer to this as the quarkyonic regime. Finally at \x d « 0.8a -1 w 850MeV 
both n q and p start to climb above the SB values. If for the moment we assume the degenerate 
quark description still holds, then for fiQ < fi < fid we have Ep « kp as expected for weakly 
interacting massless quarks, but Ep < kp for fi > fid, implying that the quark matter at very 
large densities has a large negative correction to the kinetic energy, ie. it is strongly self-bound. 

4. Order Parameters and Phase Diagram 

In order to characterise the high density system we have also examined two "order parameters" , 
one exact, the other approximate. Fig. [7] shows the superfluid order parameter (qq) as a 
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Figure 7. (qq)/fi 2 for j — > for various T. Figure 8. Polyakov line L(fi) for various T. 



function of fi, extrapolated to j = 0, for three different temperatures. For a degenerate system 
superfluidity arises through Cooper pair condensation of diquark pairs at the Fermi surface; 
hence (qq) should scale as the area of the Fermi surface kp. The plot shows (qq)/fi 2 , and 
indeed it is remarkably constant within the quarkyonic regime, although clearly T-sensitive. 
The temperature sensitivity becomes more marked once \i > fid, and at the highest temperature 
studied the order parameter vanishes everywhere, implying restoration of the "normal" phase. 

We also studied the Polyakov loop L as a function of both fi and T. Since this involves 
comparison of data with different N T , it has been necessary to renormalise L via multiplication 
by a factor Zp T [6], determined at fi = and normalised so that L{T = l/4a) = 1. The results 
are shown in Fig. [HJ the inset shows the unrenormalised data. First consider the data from 
the lowest temperature with A^ T = 24. Although in the presence of fundamental matter L is 
not an exact order parameter for global Z^ c center symmetry, its behaviour strongly suggests 
that the transition at fid is for all practical purposes identical with deconfmement, ie. the free 
energy for a static fundamental color source becomes finite once fi > fid- The quark density 
at deconfinement is 16 - 32 fm~ 3 (the uncertainty arises from the difficulty in handling lattice 
artifacts discussed in the previous section), some 30 - 60 times that of nuclear matter. 
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Figure 9. Tentative QC 2 D phase diagram. Figure 10. x q /x S q B for ja = 0.04, various T 



It's then interesting that fid defined via L falls rapidly as T rises. A pragmatic definition 
of fid(T) is the value at which L( / u,T) L(0, T^). The resulting tentative phase diagram is 
shown in Fig. [9l There are at least three distinct regions/phases: a normal hadronic phase with 



(qq) = 0, L ~ at low T and \i\ the quarkyonic region with (qq) oc fi 2 and L ~ at low T and 
intermediate and a deconfined, normal phase with (qq) = 0, L > at large T and/or large \x. 
At this stage we cannot exclude a deconfined superfluid phase at large fi and small/intermediate 
T. The "quarkyonic" nomenclature is now justified; this phase has the thermodynamic bulk 
scaling of a weakly- interacting degenerate quarks, but remains confined. It thus has two of the 
features of dense baryonic matter originally proposed on the basis of large- N c arguments [7]. 

5. Quark Number Susceptibility and Deconfinement 

We have also recently calculated the quark number susceptibility Xq = dn q /dfi [8j [3]. Whilst 
naively \ q is expected to reflect local fluctuations of the n q operator, it turns out that at large \x 
the dominant term comes from a connected "hairpin" diagram. Fig. [10] shows Xq divided by the 
continuum SB result. Once again, we note a large range over which the ratio is approximately 
constant; indeed it is compatible with one as j — > 0, though more sensitive to the quark mass 
value used for the free fermions than other bulk observables [3]. 

The most striking feature of Fig. [TOl however, is the absence of T-dependence at all but 
the highest temperature studied (141 MeV), is in stark contrast to the behaviour of L in 
Fig. [HJ Despite our intuition from the thermal transition in QCD [9], Xq cannot be regarded 
as a proxy for L once fi/T 3> 1. We conjecture that the change in behaviour of the bulk 
thermodynamic quantities n q , p and Xq observed at fid is a transition from short-ranged binary 
confining interactions to longer-ranged interactions among several quarks within the medium, 
and that this corresponds with the transition from weak to strong self-binding noted in Sec. [3j 
The T-dependent behaviour of L, in contrast, must be due to degrees of freedom close to the 
Fermi surface which can be thermally excited; a correlation between gapless excitations and 
L > has also been noted in analytical and numerical studies on small cold volumes [10j . The 
same phenomenon should inform a theory of transport properties. 

6. Summary 

QC2D offers an accessible theoretical laboratory for the study of dense baryonic matter. For low 
T at least three physical regions can be identified: the vacuum for fi < fi Q ; a confined quarkyonic 
superfluid for fiQ < (j, < fid', and a deconfined phase for fi > fid- The most recent simulations 
have only made our findings in the quarkyonic regime more robust. Not discussed here are new 
results for renormalised energy density, conformal anomaly, and chiral symmetry restoration |3] . 
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